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Abstract 

In this short note, I point out that [p, q] ^ ih, contrary to the original claims of [Born and 



Jordan, and Dirac. Rather, [p, q] is equal to something that is infinitesimally different 
from ih. While this difference is usually harmless, it does provide the solution of the 
Born- Jordan "trace paradox" of [p, q]. More recently subtleties of a very similar form 
have been found to be of fundamental importance in quantum field theory. 



When Born and Jordan ]T[ and Dirad discovered the relationship [p, q] = ih, it was 



a turning point in physics. Classically, physical quantities had always been assumed to 
commute; quantum mechanics was born when this assumption was discarded. Matrix 
mechanics reflects this non-commutativity by representing quantities such as p and q by 
matrices; wave mechanics does likewise by considering them to be operators; and Dirac's 
c-number and q-number formulation simply takes non-commutativity as the starting point. 

Born and Jordan obtained [p, q] = ih by means of arguments based on the correspondence 
principle; and Dirac obtained it by his Poisson bracket ansatz. Let us review its standard 
wave-mechanical derivation. In the g-representation, the state vector is a function ip{q), 
the operator q is simply multiplication by q, and the operator p is defined as ihd/dq. Thus 
the identity [p, q] = ih is just a scaling by the factor ih of the identity 

fo,?] = l, (1) 

where I am using the notation d q to denote d/dq. The meaning of ([!]) is made more explicit 
if we write in the implied function ip{q) on both sides: 

[d q ,q\ip(q) = ip(q)- 

It is straightforward to multiply out the left-hand side, and use the product rule to obtain 
the right-hand side: 

[dq, q]^(q) = d q qip(q) - qd q ^(q) 

= i/j(q) + qd q ip(q) - qd q ^(q) 



A problem, however, arises when we want to make the transition to the matrix mechanical 
formulation of quantum mechanics. In this formulation, the state vector is written as 



a column vector. Physically observable quantities, such as p and q, must be represented by 
Hermitian matrices. The equivalent of the g-representation of wave mechanics is obtained 
by taking the vertical position in the column vector as a linear function of the value q. Since 
q is a continuous variable, and the rows of a column vector are discrete, we must consider 
the limit of a sequence of discrete matrix representations, of ever increasing dimension, such 
that the positions in the column vector for "fill in" the domain of q more and more 
densely, so that in the limit of an infinite-dimensional matrix they form a continuum. Let 
us label the rows in such a way that the "middle one" has the index value n = 0; the rows 
above are rows n — — 1, —2, —3, . . . , and those below are n — +1, +2, +3, .... Let us then 
deem that row is to represent the origin of the q coordinate, q = 0. Then the relationship 
between q and n is of the form 

q = £n, 

where £ is some length scale, that will shrink as the dimension of the matrix representation 
is increased. (The precise mathematical form of this "shrinking rate" does not need to be 
known for our purposes). In other words, if we denote the column vector representing \%p) 
by the boldface symbol if), then we have 
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Let us now construct the matrix q that represents the operator q. Clearly, the quantity 
qip{q) is given in the matrix representation by 
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it is then clear that 
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is the matrix corresponding to the operator q, where dots indicate zero entries in the matrix. 

Constructing a matrix p to represent the operator p is a little more subtle. Clearly, it 
relies on us devising a suitable matrix operator that is equivalent to the operator d q , in 
the limit of an infinite dimensional matrix. Now, since p is postulated to be an observable 
quantity, the matrix p must be Hermitian; and since by definition p = ihd q , then it follows 
that the matrix representation d q of d q must be antz-Hermitian. Furthermore, since the 
derivative d q of any real function ip(q) must itself be real, and since the definition of d q 
cannot depend on whether the function we apply it to is real or complex, then it follows 
that d q must in full generality be a real matrix. Taken together, these two considerations 
already tell us that d q must be a real, antisymmetric matrix. To find its exact form, let us 
consider the meaning of the derivative d q from first principles: for a function ip(q), 
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(3) 



where e and e' are real numbers greater than zero. Now, in the matrix representation, for a 
finite dimension, we do not have positions that are infinitesimally close to a given q n = n£; 
rather, the closest we can get are the two points q n+ \ = (n + 1)£ and = (n — 1)£. 
However, in the limit of an infinite-dimensional matrix representation, these two points will 
shrink around the point q n in the way we desire. Moreover, we already know that we cannot 
use the point q n itself in the definition of d q , since the matrix d q must be antisymmetric, 
which means that the diagonal elements must vanish. The best that we can do is therefore 



(d q if>) n 
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which is equivalent to ([3D, with e = e f , in the limit of an infinite-dimensional matrix. This 
then imples that 
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as can be verified directly by multiplying ((^) by (0). In other words, the matrix p is given 
by 
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which is clearly Hermitian, as required. 

We can now turn immediately to the issue raised by the title of this note, by computing 
[p, q] in the matrix representation above — namely, by computing the matrix commutator 
[p, q] . By multiplying out the matrices, it is easily seen that 
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we therefore find that 
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Here is the subtlety. The problem is that the matrix (|7]) is not equal to ih times the unit 
matrix 1, 



ih l = — 

2 
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Rather, the matrix (0) is effectively obtained by taking each diagonal element of 2 and 
"splitting it" between the off-diagonals above and below. Thus we have proved the relation 



[p, q]^ihl 



(9) 
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in the matrix representation of quantum mechanics, and hence in full generality 



\p,q]^ih, (10) 

as I have claimed in the title of this note. 

The result (|10"D might be somewhat disturbing. However, in almost all cases, it is of 
academic interest only. The reason is that using the matrix (§) rather than fl7|), in any 
practical calculation, corresponds to the replacement 

e->o 2 

which is arguably harmless for any reasonable ip(q). In fact, we can obtain exactly the 
same result [p, q] ^ ih using the wave mechanical representation, if we treat the operation of 
differentiation more carefully, rather than by simply using the product rule. If we write down 
the wave-mechanical equivalent of the Hermitian (i.e., symmetrical) definition (f|), namely, 
the symmetrical version of 



d q 4{q') 
then we find that 

[d g , q}^(q) = d q qtp(q) - qd q il)(q) 



= lim ^(g + g) -ip{g-e) 
g'=g ~ e^o 2e 



~ ™\ 2e q 2e 

if>(q + e) + ip(q - e) 



= lim 

in agreement with the matrix mechanical result. 

It might seem that claiming that [p, q] ^ ih is a pedantry. After all, when would shifting 
the argument q by an infinitesimal amount, or shifting by one row or column in an infinite- 
dimensional matrix representation, make any difference? There is at least one situation that 
I am aware of in which this change does make a difference: whenever the trace of the matrix 
is taken. For example, the Born and Jordan's |1[ well-known "trace paradox" of [p, q] points 
out the following: since 

Ti(AB) = Ti(BA) 

for any finite matrices A and B, then in the finite-dimensional case we must have 

Tr[p, q ] = Ti(pq-q P ) = 0. 
But if [p, q] = ih were to hold true, then we would need to have 

Tr[p,q] = zftTr(l) = ih D, 

where D is the dimension of the matrix representation, which, rather than vanishing, ap- 
proaches infinity in the infinite-dimensional limit! I emphasise that this is a fallacy; it is 
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the matrix (0) that must be used, not the identity matrix (H). And of course the matrix 
(0) is identically traceless; hence, the Born- Jordan "trace paradox" of [p,q] is due to the 
incorrect assumption that [p, q] = ih, whereas at the level of individual rows and columns of 
the matrix representation it fails. 

It might be claimed that this simply shows that one cannot take infinite-dimensional 
matrix mechanics to be the infinite-dimension limit of finite-dimensional matrix mechanics. 
But then what would this "matrix mechanics" have to do with matrices as we know them? 
Moreover, it is generally believed that the correct way of dealing with infinities, or infinites- 
imals, in physical problems is to take them to be the limit of finite quantities. Surely, then, 
it is better to modify the postulate of [p, q] = ih by an infinitesimal amount, rather than 
remove all chance of using a well-defined limiting procedure? 

Furthermore, this ability of the trace — to be able to yield an answer that is either zero 
or infinite, depending on how carelessly one defines one's matrix quantities — turns out to be 
more important to real-world calculations than one might naively think. In quantum field 
theory, the effect of effectively "including the diagonal terms" in the time-ordered product 
operation, when in fact they should not be included, leads to a drastic and fundamental 
change in the predictions of calculations involving loop diagrams. This has has been pointed 
out several times in the past two decades, but has not gained much attention; we shall be 
providing a full and explicit description of these developments shortly [[| . 
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